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ABSTRACT 



. Aims. We investigate constraints on additional planets orbiting the distant M-dwarf star OGLE 2005-BLG-390L, around which 

photometric microlensing data has revealed the existence of the sub-Neptune-mass planet OGLE 2005-BLG-390Lb. We specifically 
aim to study potential Jovian companions and compare our findings with predictions from core-accretion and disc-instability models 
of planet formation. We also obtain an estimate of the detection probability for sub-Neptune mass planets similar to OGLE 2005- 
BLG-390Lb using a simplified simulation of a microlensing experiment. 

Methods. We compute the efficiency of our photometric data for detecting additional planets around OGLE 2005-BLG-390L, as 
a function of the microlensing model parameters and convert it into a function of the orbital axis and planet mass by means of an 
adopted model of the Milky Way. 

Results. We find that more than 50% of potential planets with a mass in excess of 1 Mj between 1.1 and 2.3 AU around OGLE 2005- 
BLG-390L would have revealed their existence, whereas for gas giants above 3 Mj in orbits between 1.5 and 2.2 AU, the detection 
efficiency reaches 70%; however, no such companion was observed. Our photometric microlensing data therefore do not contradict 
the existence of gas giant planets at any separation orbiting OGLE 2005-BLG-390L. Furthermore we find a detection probability for 
an OGLE 2005-BLG-390Lb-like planet of around 2 — 5%. In agreement with current planet formation theories, this quantitatively 
supports the prediction that sub-Neptune mass planets are common around low-mass stars. 

Key words, techniques: microlensing, exoplanets: individual: OGLE 2005-BLG-390Lb, stars: M-dwarfs, K-dwarfs 



1. Introduction 
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After more than a decade since the disco very of the first 
extra-solar planet orbiting a solar-type star dMavor & Ouelozl 
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1 1995b the Extra-solar Planet Encyclopedia^ lists 249 entries, in- 
cluding an increasing number (26) of multiple planetary sys- 
tems. The accessible mass regime of extra-solar planets ex- 
tends since 2005 below the 10 M y regime, with the discov- 
eries of Gliese 876 d (~ 7.5 M e . iRivera et alj|2005l) and the 
pair of pla nets around Glie se 581 with minimum masses of 5 
and 8 M m dUdrv et alj2007l) using the radial-velocity technique, 
as well as OG LE 2005-BLG-390Lb detected by microlens- 
ing (~ 5 M®, iBeaulieu et al.l I2006T) . While recent improve- 
ments in radial velocity sensitivity have ena bled the discovery 
of Neptune-mass p lanets in Venus-like orbits (lLovis et al ] |2006t 
lAlibert et al.l2006l) . microlensing is the only method that can de- 
tect such sub-Neptune-mass planets in orbits beyond 1 AU. This 
sensitivity to sub-Neptune-mass planets at separations of a few 
AU is important for testing the core accretion theory of planet 
formation because this theory predicts that the dominant plan- 
ets in any planetary system should form in the vicinity of the 
"snow line',' which is located at a few AU dKennedv et alj|2006l 
llda & Linl 120041: ILaughlin et"ai1l2004l) . Microlensing results al- 
low this theory to be tested without confronting the additional 
uncertainties of planetary migration. 

Despite these impressive successes, the observational picture 
of the planet abundance and their mass function is far from com- 
plete. This is partly caused by biases introduced by the detec- 
tion techniques. The currently dominant radial-velocity method, 
despite ever improving sensitivity and temporal baselines, still 
favors massive planets in close-in orbits around solar type stars 
within about 100 pc from the Sun. The transit method also favors 
the detection of close-in giant planets. 

The microle n sing technique, initially proposed by 
iMao & Paczvnskil d!991l) and whose prospects h ave been 
quantified first in more detail bv lGould & Loebld 19921). has now 
led to four reported detect i ons dBond et al. 2004 ; lUdalski et all 
|200H IBeaulieu et alj|200d iGould et al.ll2006l)~ It has proved its 
capability to provide access to a new window for exoplanets, 
with masses down to an Earth-mass for ground-based searches 
and orbits between 1 — 10 AU mainly around host stars less 
massive than the Sun at several kpc distance. This is of special 
interest since most of stars in our Galaxy have masses less than 
1 Mq, so the planets of such stars, if common, may constitute 
the majority of Galactic exoplanetary systems. While the 
hunt for planets around low-mass stars still lags a bit behind 
the search for planetary companions to solar-type stars, the 
most recent microlen sing discoveries of sub-Neptune mass (or 
super-Earth) planets dBeaulieu et all 120061: IGould et alj 120061) 
indicate that such planets indeed are common, probably more 
common than any other class of exoplanets yet discovered. The 
rarity of cool Jovian companions to sub-so lar mass stars (from 
1 — 10 AU) seen in the microlensing data dAlbrow et al]|200U 



iGaudi et alJl2002l:lTsapras et al.l2003tlSnodgrass et alj|2004 on 
the other hand was recently complemented by radial velocity 
searches, deriving an upper limit of < 1% for the fraction of 
close -in (within 1 AU) Jovian M-dwarf planets (Endletal. 
120061) . 

These new observational constraints seem to be in line with 
predictions o f planet formation the ories. For example sim ula- 
tions done by llda & Linl d2005l) and ILaughlin et al] d2004l) sug- 
gest that the formation of Jovian gas giants around M-dwarfs 
is inhibited, while planets less massive than Neptune can easily 
form. This is supported by recent microlensing results as one of 
the gas-giant planets disc overed by microlens ing was recently 
found to orbit a K-dwarf dBennett et al.ll2006l) . while the least 



massive planets discovered so far (OGLE 2005-BLG-390Lb, 
Gliese 476 b, Gliese 581 c and d) are all orbiting M-dwarf stars . 
Disc instability formation models as advocated by iBossI d2006l) 
also are capable of explaining the preference of forming sub- 
Neptune mass planets rather than giant planets around low-mass 
stars . 

In the study presented here we examine what constraints on 
a hypothetical Jovian planetary companion to the sub-Neptune 
mass plan et detected in the mi crolensing event OGLE 2005- 
BLG-390 dBeaulieu et al.ll2006b can be derived from the pho- 
tometric light curve data. This system contains a 3 — 10 M® 
planet in orbit with semi major axis a = 2 — 5 AU around an 
M-dwarf. We also estimate the probability of having detected an 
OGLE 2005-BLG-390Lb like planet in an idealized microlens- 
ing experimen t and d iscu ss in more details t he claims made in 
IBeaulieu etaf] (120061) and lGould et al.ld2006l) that sub-Neptunes 
are common companions to low-mass stars. 



2. Basics of microlensing 

During a microlensing event, the light arising from a background 
star (the source star) is bent due to the gravitational field of 
an intervening planetary system (the lens) passing close to the 
observer line-of-sight. This results in a characteristic transient 
magnification of the source leading to a variation in the received 
flux, which constitutes the microlensing event light curve. 

With M denoting the total lens mass and Dg, £>l the dis- 
tances of source and l ens from the ob server respectively, the an- 
gular Einstein radius dEinsteinll 19361) 



4GM fD s -D L 



(1) 
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defines the natural scale of gravitational microlensing. It equals 
the angular radius of the ring-shaped image of the source star 
which would occur for a perfect alignment between the observer, 
lens and source. With sources typically located in the Galactic 
Bulge and a lens lying at several kpc from us in the direction 
of the source, its linear scale at the lens distance Re = -Dl #e 
corresponds to ~ 1 — 10 AU, a range well suited for extrasolar 
planet hunting. 

The transient brightening of the source (called magnifica- 
tion) due to a single lens can be fairly easily modeled and is 
described by an impact parameter uq, with uo #e being the min- 
imal angular distance between lens and source at time to, as well 
as a time-scale ig which is the duration for the source to move 
by #e on the sky. The total observed flux F(t) is then expressed 
as the sum of the magnified source flux A(t) x F$, with A(t) 
being the magnification factor, plus the so-called blend flux Fb 
stemming from unresolved sources within the aperture, so that 
F(t) = A(t) Fs + Fb- Furthermore, one must take into account 
a pair of flux parameters (Fg, Fb) per observing telescope. 
A lens with a planet is usually modeled by a binary lens with 
an extreme mass ratio (q < 10 -2 ). The most striking char- 
acteristic of binary lenses when compared to single lenses is 
the occurrence of extended caustics, closed lines defined in the 
source plane where the magnification becomes infinite for point 
sources. However in practice the magnification remains finite 
since real sources have a finite extent. Nevertheless caustics 
mark regions of large magnification gradients. In the case of 
planetary lensing the area enclosed by these caustics (defining 
a kind of cross section) compared to the area covered by the 
Einstein ring disc is only a few percent at most, but if the source 



D. Kubas, A. Cassan et al. : Multiplicity of OGLE 2005-BLG-390L 



magnification [mag] 



nonlinear superposition 



0.3 h 
0.2 
0.1 - 
0.0 
-0.1 



Jovian planetary caustic 



central caustic 




OGLE-20u5-BLG-390Lb 
planetary caustic 



-0.2 



0.0 



0.2 



0.4 

V 6 E 





Fig. 1. The left panel of the figure displays a magnification map of a triple lens system configuration, consisting of a host star at 
the origin and two planets. The magnification is shown as function of the position angle [6\, #2) on the sky. The darker the region 
the higher the magnification. The caustic contours have been over-plotted (black solid lines) for clarity. The small planetary caustic 
induced by OGLE 2005-BLG-390Lb is shown together with the best-fitting source track (iBeaulieu et al.ll2006l pink solid line). 
The larger planetary caustic on the upper left stems from a hypothetical additional Jovian companion, in this particular case with 
q = 2.0 x 10~ 3 and d = 1.2. Note that the actual positions of the planets do not coincide with the position of their planetary 
caustics. They lie on the dotted axes but outside the shown field of view, which scale was chosen to maximize the visibility of the 
caustics. Finally the right panel shows one of the rare configurations where the caustics of the two planets are close to the same 
position merging in a non-linear way, showing the limit of the linear superposition approximation. 



trajectory passes sufficiently close to or even traverses a caus- 
tic it can imprint a high signal-to-noise light curve signature, 
revealing the presence of the planet. Unfortunately this signal 
is only short-lived (ranging from hours to days for Earths to 
Jupiters) with respect to the complete event time-scales of weeks 
to months. This is the big observational challenge for teams like 
PLANET/RoboNefl OGLE3, MOA0 and MicroFUI\0 to moni- 
tor microlensing events with a high sampling rate and not to miss 
these planetary 'anomalies' in light curves which for most of the 
time are indistinguishable from a single-lens event. 

The challenge on the modeling side is that the magnification 
of a binary lens cannot be expressed in a closed analytical form 
and its high-dimensional parameter space has an intricate x 2 - 
surface on which the problem of parameter optimization is far 
from trivial. In addition to the single-lens model parameters the 
following ones are required: the planet-to-star mass ratio q, the 
lens separation d, with d x 9e being the instantaneous angular 
distance of the planet from its host star, the impact angle a be- 
tween the source trajectory and the binary lens axis. If the source 
is resolved by the lens, the angular source size 0* and parameters 
describing the source surface profile also need to be considered. 

3. Detection efficiency of additional planets 
around OGLE 2005-BLG-390L 

iGaudi & Sackettl d2000l) have presented an algorithm well-suited 
to compute the "efficiency" with which a binary companion to 
the lens can be revealed in an observed microlensing event. This 

2 http://planet.iap.fr 
http://www.astro.livjm.ac.uk/RoboNet/ 

3 http://www.astrouw.edu.pl/~ogle 

4 http://www.phys.canterbury.ac.nz/moa 

3 http://www-astronomy.mps. ohio-state.edu/~microfun 



algorithm applies for microlensing events where no clear de- 
viation from a single point-mass lens can be seen in the light 
curve, i.e. events for which there is no evidence f or a compan- 
ion to the lens. Efficiency ca lculations as done in lAlbrow et alj 
d2000b or lGaudi et aT] i]2002l) indeed were aimed to give limits 
on companions to lenses where no planets were found. Here we 
extend this method to compute the detection efficiency of multi- 
ple planetary systems, more precisely to put confidence limits on 
the presence of an hypothetical second planet in a system where 
one planet has already been detected. For similar reasons to the 
original method, the algorithm holds as long as there is no fur- 
ther clear signal in the light curve other than the already detected 
planet. We first outline the method before applying it specifically 
to OGLE 2005-BLG-390L. 

As a natural exte nsion of the method proposed by 
IGaudi & Sackettl (|2000|) . we use a triple lens configuration to 
describe the discovered planet plus and additional companion, 
for which we want to compute detection efficiencies (i.e. calcu- 
lating confidence limits on the latter's existence). Such a geo- 
metrical configuration is shown in the left panel of Fig. 1, where 
the parent star of the system is located at the origin and known 
OGLE 2005-BLG-390Lb is situated on the horizontal axis on 
the right. A hypothetical additional planet with mass ratio q is lo- 
cated at a distance d from its parent star (d = 1.2, q = 2 x 10~ 3 
and <j> = 60° in figure) with its planet-star axis subtending an 
angle (f> with the .T-axis. Note that with a planetary signal al- 
ready detected, the trajectory of the source is known and fixed. 
We then compute detection efficiency on the second planet using 
the algorithm described below: 

1. Compute the x^ e{ of the model that best fits the already dis- 
covered planet, 

2. Choose a given triple lens configuration with parameters 
(d, q, </>), where <fi <E [0, 27r], and fit for the remaining 
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parameters to compute the \ 2 difference A\ 2 (d, q, <fi) = 
X 2 (d,q,(f>) - X?ef- 

3. Repeat step (2) until a dense sampling of <fi is obtained for 
the probed lens configurations, 

4. Repeat from step (3) until all the chosen lens configurations 
are probed. 

The detection efficiency of an additional planet in the system 
is then given by the fraction of angles 4> that produce a signifi- 
cant deviation in the light curve. By significant deviation, we 
mean an excess of x 2 of the probed lens model by an amount of 
Xthrcsh restive to the single planetary model. More formally the 
detection efficiency is given by the formula: 

<d, q) = j H {A x 2 (d, q, 0) - xLcsh) (2) 



H being the Heaviside step function. 

Apart from the computationally challenging sheer size of the 
parameter space of multiple lens systems another notable diffi- 
culty arises when the spatial extent of the source star has to be 
considered. Rather than integrating over the source surface, we 
make use of magnification maps that are much better designed 
for our purpose. A magnification map is obtained by shooting 
light ra ys from the observe r through the lens and towards the 
source dWambsgansslll997h ; assuming a constant spatial ray 
density in the plane containing the lens, its density in the plane 
where resides the source then maps the magnification as a func- 
tion of its position. Such maps can easily include the extended 
source effects: The ray-shot maps just need to be convolved us- 
ing the ri ght size and brigh tness profile of the chosen source 
star (e.g. iKubas et al]|2005l) . A speed optimiz ed version of the 
origin al ray-shooting algorithm is presented in lRattenburv et all 
d2002l) . Fig. 1 is an example of such a triple lens magnification 
map. 

In general the source size in physical or Einstein units is not 
known and estimates rely on statistical Galactic models. The un- 
certainty of those parameters directly affec ts the uncertaintie s 
of the computed detection efficiencies (cf. iGaudi et"aTll2002[) . 
The situation is much more favorable when, as in the case of 
OGLE 2005-BLG-390Lb, extended source effects are detected 
and strong constraints can be put on the source properties. Then 
it is straightforward to convolve a library of maps with the ap- 
propiate size and brightness profile of the source which as a wel- 
comed byproduct reduces the to be explored parameter space. 

Nevertheless this approach is still resource intensive since 
it requires computing a large number of magnification maps for 
the different probed characteristics (d, q, </>) of the second planet. 
However, in the case of OGLE 2005-BLG-390 we can simplify 
the problem with a fair approximation as follows. In fact, it has 
been shown that the planetary caustics in a multiple planetary 
system are independent to firs t order, as lo ng as the projected 
planet positions do not overlap dBozzal 1999b . However when the 
planetary caustics of the two planets are too close, non-linearities 
appear and the superposition principle is not valid anymore. This 
effect can be seen in the right panel of Fig. 1, where the two 
caustics are superposed in a non-linear way. Here the more mas- 
sive hypothetical planet has been located only 0.05 fractions of 
i?E above OGLE 2005-BLG-390Lb causing the small caustic 
to flip its orientation by about 90 degrees and splitting two of its 
cusps. Hence, the superposition principle of the caustics is valid 
for most values of the angle cf> and position d of the additional 
planet, except when the caustics coincide. We make the choice 
here to compute only binary lens maps and use the superposition 



principle. By doing so, we save around two orders of magnitude 
of computational time to create the map library (each magnifica- 
tion map requires few hours of computation with a usual CPU). 

To deal with the rare configurations when the non-linearities 
appear, one choice would be to subtract the best-fit binary-lens 
model from the original OGLE 2005-BLG-390Lb light curve 
and then add the residuals to form a new light curve based on 
the underlying single-lens model. This new "single-lens" light 
curve then serves as a referen ce to which a series of b inary-lens 
models are fitted according to Gaudi & Sackett (2000). However 
we prefer here to adopt a more conservative choice, by cutting 
out the data from the region where the signal of OGLE 2005- 
BLG-390Lb resides (i.e MHJD between 3592.50 and 3593.11) 
prior to the computations. While with both options one cannot 
remove the intrinsic non-linearity at the common position of the 
caustics, the uncertainties they introduce are small and in any 
case of similar order of magnitude. 

We have computed a (d, q) magnification map library of 
lens configurations spanning nine different mass ratios (q = 
5.6 x 10~ 4 — 10~ 2 ) and covering sixteen lens separations (d = 
0.1 — 3.3), convolved wi th the source brightness profile adopted 
by iBeaulieu et al.l ([2006). The choice of the grid is motivated 
by the goal to cover the theo retical range of Jo vian type planets 
probed in the simulations of llda & Linl d2005l) . and also to en- 
sure a good coverage of the so-called lensing zone, the range of 
lens separations, where microlensing is most sensitive to plan- 
ets. For each (d, q) configuration we have computed 200 models 
with 4> k = 2fc7r/200, i.e. in total 200 x 16 x 9 = 28800 models. 

At each (d, g)-grid point and for each angle (f> 6 [0, 2tt] 
we then compute the least-square measure x 2 (d, q, <t>) optimized 
over the remaining free parameters: minimum impact parameter 
uq, time of maximum magnification to, Einstein radius cross- 
ing time tE, as well as the source fluxes F^ 1 and blending fluxes 
F-q for each of the 6 different telescope data sets used. These 
optimized parameters are 'free' in the sense that we allow them 
to vary within th e error bars of the best-fitting values derived 
in IBeaulieu et al.l (120061) . This flexibility minimizes numerical 
noise in the calculations without violating the constraints given 
by the best model. We checked that the results are consistent 
with fixing the parameters to their best-fitting values and that us- 
ing non-bounded parameters yields unphysical re sults. To carry 
out th e optimization we use a genetics algorithm dCharbonneaul 
fl995l) . which naturally provides the capability to bound parame- 
ters and has been shown to explore the intricate parameter space 
of binary lenses mo re efficiently than classica l gradient opti- 
mization techniques (Kubas 20051: ICassanll2005l) . 

A given choice of parameters (d, q, <fi) is considered to pro- 
duce a significant deviation if Ax 2 = X 2 (d, q, 4>) ~ X?ef ex ~ 
ceeds a threshold value Xthrcsh- We set Xthrcsh = 60, which we 
find robust enough to avoid false detections arising from statisti- 
cal fluctuations or unrecognized low-level systemati cs and to be 
consistent with earlier detection efficiency studies ( IGaudi et~aT1 
l2002l:lYoo etaill2004l) . 

Fig. [2] shows contours of the detection efficiency e(d, q) as a 
function of the dimensionless model parameters d and q, where 
d x 6-e is the angular separation of the planet from its host star 
while q is the planet-to-star mass ratio. Calculations using an al- 
gorithm based on a method described in iRhie et al.l d2000l) were 
also used to compute detection efficiencies, with comparable re- 
sults. 

The obtained detection efficiency constraints are admittedly 
rather weak. However this is not a flaw of the data or analysis but 
mainly due to the intrinsic nature of a low magnification event. 
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iBeaulieu et all (l2006l) . namely D s ~ 8.0 kpc, £> L ^ 6.6 kpc 
and M ~ 0.22 M a , we have: 



lens seporotion d (0 E ) 

Fig. 2. Detection efficiency e(d,q) as a function of angular 
planet-to-star separation d x E and planet-to-star mass ratio g. 
These dimensionless efficiencies present the raw outcome of our 
computation prior to the convolution given in Appendix A where 
we derive the physical detection efficiency values as functions of 
planet mass and orbital radius. 



In fact, iGaudi & Sacked d2000h have shown that the detection 
efficiency of a microlensing event is strongly dependent on the 
minimum parameter uq of the underlying single lens model: the 
smaller uq, the higher the peak magnification and the higher the 
detection efficiency. This can be understood as follows. 

As depicted in Fig. 1, in a planetary binary lens scenario one 
has to differentiate between the central caustic, always located 
close to the primary lens (star), and the planetary caustics, the 
location of which strongly depends on the projected star-planet 
separation d. While as stated earlier the planetary caustics are ba- 
sically mutually independent from each other, the central caus- 
tic is affected by any companion to the primary lens and thus 
also especially sensitive to multiple planet systems. To probe 
the central caustic sufficiently small impact parameters up, re- 
spectiv ely high magnifications, are required fGriest & Safizadehl 
(119981) showed that for uq — » the planet detection efficiency 
goes to one, unless finite-source effects prevent the detection 
of the least massive planets around the highest peaks for larger 
sources. The drawback of high magnification central caustic 
events however is that they are much harder to model and some- 
times plagued by a close-wide binary ambiguity dDomini kll999t 
iKubas et al.ll2005h . In the OGLE 2005-BLG-390 event the im- 
pact parameter is too large to explore the central caustic, which 
explains the calculated restricted sensitivity to additional plan- 
ets. 



4. Limits on Jovian companions to 
OGLE 2005-BLG-390L 

While E m (d,q) is the most straightforward constraint we can 
get from modeling, one would like to infer statements on the 
underlying physical parameters, the planet's orbital axis a and 
its mass m. The planet mass and the mass ratio are linked by 
m = q X M, with M being the host star mass. The phys- 
ical (instantaneous) lens separation projection r is given by 
r = d x Dl x 6>e- Adopting the median values determined by 



m = qx M ~qx 0.22 M 
r = dx D L 6» E ^dx 1.35 AU. 



(3) 



However, without a proper measurement of both the source 
size and the para llax in Einstein radii (eg. iGhosh et al.l 120041: 
Uiang etal.H"2004l) . the mass M of the host star OGLE 2005- 
BLG-390L and the Einstein radius are distributed in probabil- 
ity over a finite range, where both distributions are correlated, 
and moreover, the orbital axis follows from r by a stochastic 
orbital de-projection. Therefore, one needs to take these distri- 
butions into account for expressing the detection efficiency as 
function of (a , to). Assuming c ircular orbits and the stellar mass 
function from lChabriei ( 20031). we use the implementation of a 
Galactic model by iDominikl J2006I) in order to derive detection 
efficiency values e(a, m). The details of this procedure are given 
in Appendix A. 

The resulting detection efficiency diagram is presented in 
Fig. [3] which shows a contour plot of e(d, q) together with 
OGLE 2005-BLG-390Lb (red error bar cross). It basically 
tells us about confidence limits on additional companions to 
OGLE 2005-BLG-390Lb: values of e(d, q) close to unity rule 
out the possibility of such an additional planet at (a, to), while 
values close to zero mean that no conclusion can be drawn about 
the multiplicity of OGLE 2005-BLG-390L. 

The strongest constraints from the diagram are placed in the 
Jovian-mass regime, where e(a, m) > 50 %. We find that plan- 
ets more massive than 1 Mj in the orbital range 1.1 — 2.3 AU 
have a detection efficiency in excess of 50 %, while planets 
above 3 Mj in orbits between 1.5 and 2.2 AU would have re- 
vealed their existence with a probability of more than 70 %. 
Unfortunately, our data cannot tell us about possible Pegasid 
planets (giant planets orbiting at fractions of 1 AU), which could 
exist in OGLE 2005-BLG-390L. 

On the same Fig . [3] we have plot ted the results of Monte- 
Carlo simulations bv llda & Linl (12005). producing the final evo- 
lution stage of a seed of 20000 planetary embryos, uniformly 
distributed in logo (from 0.1 to 100 AU), around a host star 
with M = 0.2 Mq, similar to OGLE 2005-BLG-390L. In 
principle, such simulations are valid when considering a single 
planet, whereas we consider the case of t wo planets . We th ere- 
fore choose among the available models of llda & Linl d2005l) one 
for which the planet migration process is very inefficient and so 
practically assume that the planets were formed quasi-wi situ. 
This allows us to use the efficiency diagram to compare our ob- 
servation to this planet formation model. I n the chosen plan et 
formation scenario - model from Fig. 9b of llda & Linl d2005b - 
migration is strongly suppressed, the cores have more time for 
accretion, more gas giants can form and planets (blue points on 
Fig- Ell stay close to their orbital birth places. Two main predic- 
tions can be read out from their simulation. Firstly, sub-Neptune- 
M-dwarf planets should vastly outnumber gas-giant-M-dwarf 
ones. In the case of an 0.2 Mq M-dwarf the formation of a sub- 
Neptune planet between 1 and 14 M® in the orbital range of 
1-10 AU is about ~ 1200 times more likely than forming Jovian 
planets with 0.5 — 1 Mj in the same orbital range. Secondly if gi- 
ant planets form they are unlikely to become more massive than 
~ 1 Mj. 

In the framework of disc instability planet formation theory, 
the preference for sub-Neptune planets can also be explained by 
photo evaporation of th e gas envelo pes of giant protoplanets, as 
recently pointed out bv lBossI (120061) . However in this case there 



6 



D. Kubas, A. Cassan et al: Multiplicity of OGLE 2005-BLG-390L 



are no obvious reasons for postulating a gap in the planetary 
mass function between Jupiter and Neptune masses. 
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Fig. 3. Detection efficiency i p (a, m) for additional planets or- 
biting OGLE 2005-BLG-390L as function of orbital separation 
a and planet mass m, where contours at 5, 20, 50 and 70 % are 
shown. The cross marks the median values for the properties 
of OGLE 2005-BLG-390Lb along with 68.3% confidence inter- 
vals and the dashed horizontal lines mark the masses of Jupiter 
(J), Saturn (S), Neptune (N) and Uranus (U) for comparison. 
The blue dots represent the predicted final distribution of a seed 
of 20000 planetary cores around an M-dwarf of 0.2M Q result- 
ing from a core-accre tion model assum ing inefficient migration 
(taken from Fig. 9b of llda & Lid d2005l) ). 



We note that also around the other known low-mass planet 
found by microlensing OGLE 2005-BLG-169Lb, which host- 
star is most likely a K-dwarf, there were no traces found in the 
data of additional Jovian c ompanions in the orbital range be- 
tween about 0.5 and 14 AU dGould et al.ll2006l) . 

The current sample of known planets around M- and K- 
dwarfs is still too small to distinguish between the alternative 
scenarios of planet formation and migration. However, just a few 
more detections are likely to add valuable further information. 

5. On the detection of 

OGLE 2005-BLG-390Lb-like planets 

iBeaulieu etaT] d2006[) and iGould et alj d2006l) have stated that 
their respective detections of the low-mass planets OGLE 2005- 
BLG-390Lb and OGLE 2005-BLG-169Lb imply a large abun- 
dance of their siblings. The sensitivity of a microlensing light 
curve to the presence of planets depends strongly on, besides 
data sampling and quality, the event magnification and to some 
extent on the event time-scale and source size as well. Increasing 
the angular size of the source has, for instance, the dual effect of 
increasing the duration of the planetary signal while, at the same 
time, decreasing its amplitude. The latter effect ultimately re- 
stricts the capabilities for detection of planets to ~ 1 M s for 
giant s ources and to ~ 0.1 M m for turnoff stars in the Galactic 
bulge dBennett&Rhidll996[) . 

An accurate determination of planet abundances therefore 
requires an analysis of the detection efficiency of a representa- 
tive sample of the whole experimental data set and a comparison 



with the actual detections. For the PL ANET campaign, th is will 
be worked out in an upcoming study (C assan et al. 2008). Here, 
we run simulations to derive estimates of the probability for de- 
tection of sub-Neptune or Jovian planets based on simplified as- 
sumptions. 

We choose to pick impact parameters uq from a uniform dis- 
tribution, which is close to what is actually realized in nature ; 
however, only those events with a magnification above a char- 
acteristic threshold will be observable. Neglecting blending, we 
sample uq E [0, 1] in 100 equally spaced steps, corresponding 
to magnifications A > 1.34. While this means that we practi- 
cally cover the whole range of magnifications in the alerted and 
followed-up microlensing events, we neglect the preference for 
higher magnifications in the follow-up campaigns. 

In Sec.[3]we made use of a A\ 2 -based criterion for ruling out 
additional planets in the OGLE 2005-BLG-390L system. While 
in principle such a criterion works for detections as well we note 
that in practice one often finds that even if this criterion is ful- 
filled, a lack in coverage and/or precision makes it impossible to 
characterize the planetary model unambiguously. Unlike a rejec- 
tion a convincing detection has to meet stronger criteria (as, for 
instance, not showing systematic trend s and having a sufficient 
number of data points dVerma ak 2000)), which are not reflected 
in a criterion solely based on x 2 - In other words we want to avoid 
cases in which planets are detectable but not characterizable, re- 
spectively cannot be distinguished from non-planetary solutions. 

To minimize the effect of this caveat on our simulation we 
therefore adopt a more demanding criterion to calculate the ef- 
ficiency to discover OGLE 2005-BLG-390Lb-like planets. For 
a "discovery", we demand that the planetary signal amplitude 
exceeds 2% for at least 15 measurements, assuming an hourly 
sampling and a photometric accuracy of 1%. Neglecting non- 
white noise this criterion then translates into Ax 2 = 60, i.e. is 
consistent with the detection threshold chosen in Sec. [3] This 
criterion ensures that the planetary light curve signature is also 
well sampled, which is essential for a proper characterization of 
the lens system. For comparison we also compute the efficien- 
cies for the following set of criteria: The signal amplitude should 
exceed 2% for at least 10 and 25 measurements. 

Since the detection efficiency reaches a maximum in the so- 
called "lensing zone", corresponding to a star-planet separation 
range of 0.62 < d < 1.62, we can restrict our sampling to the 
time range [to — tE,io + <e]> as most microlensing campaigns 
do, which ensures being sensitive to planets at these separations. 
To avoid border effects, we however draw the simulated light 
curves from magnification maps spanning to ± 1.2 1^. For the 
source star, we adopt an angular size (9* = 9.6 x 10~ 3 6>e (which 
corresponds to 8* = 5.25 //as or a physical radius i?* ~ 10 Rq 
in the case of OGLE 2005-BLG-390Lb). 

As in Sect. 3, we use the ray-shooting technique for calcu- 
lating magnification maps. While we adopt fixed mass ratios of 
q = 7.6 x 10~ 5 , which is the value for OGLE 2005-BLG- 
390Lb, and q = 4.3 x 10~ 3 , which corresponds to M = 1 Mj if 
the mass of OGLE 2005-BLG-390Lb is 5.5 M B , we use a grid 
for the lens separation d = 0.1 — 5.0. Again as in Sec. [3] we 
compute the fraction of detections made over this grid assuming 
that every lens has such a planet by averaging over the impact 
angle and impact parameter. 

In Fig.|4]the derived efficiencies are shown as a function of 
the projected separation d, for the selected mass ratios and time 
scale, < E = 11.0 as measured for OGLE 2005-BLG-390. For 
criterion C2 (more than 15 deviating points) at the separation 
parameter for OGLE 2005-BLG-390Lb, d = 1.61, the detec- 
tion efficiency for an OGLE 2005-BLG-390Lb-like mass ratio 
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Fig. 4. Efficiency of detecting a super-Earth planet with the 
mass ratio of OGLE 2005-BLG-390Lb (? = 7.6 x 1(T 5 , lower 
set of curves), and for a giant planet with a Jupiter/Sun mass ra- 
tio (q = 4.3 x 10 -3 , upper set of curves) as a function of the 
projected planet-star separation d. The bold black curves have 
been computed assuming an hourly sampled light curve with 
1% photometric accuracy, and demanding that the planetary sig- 
nature amplitude exceeds 2% for at least 10 measurements (dot- 
ted line, indexed as Ci), 15 and 25 measurements (solid line, 
C2 and dashed line, C3). For comparison efficiencies for the 
same sampling rate and accuracy but using a pure A\ 2 threshold 
are plotted as well, with Ax 2 = 60, 220, 400 (blue thin lines). 
The considered impact parameters are drawn uniformly from 
uq E [0, 1]. Based on the parameters found for the lens sce- 
nario of OGLE 2005-BLG-390, a time scale of t E = 11.0 days 
was adopted, as well as a source size of ~ 10 Rq. The vertical 
line marks the actual separation of OGLE 2005-BLG-390Lb. 



events followed up by planet, the simulations yield that Jovians 
are about 10 times easier to detect that sub-Neptunes. 

However, there are some notable differences between the 
assumptions for our simulation and the real campaigns. First, 
even in dense telescope networks, weather and technical losses 
are unavoidable and a continous data coverage can not always 
be achieved. Secondly the peak magnification of OGLE 2005- 
BLG-390Lb (A a ~ 3, corresponding to uq = 0.36) is one of 
the lowest among the events monitored by PLANET/RoboNet. 
Follow-up teams preferentially monitor events with larger mag- 
nifications, thereby preferring smaller uq. As stated at the end of 
Sec. 3, the detection efficiency usually increases with the magni- 
fication and therefore the calculated efficiency for uq uniformly 
chosen in [0, 1] underestimates the chance to detect less-massive 
planets in high magnification events. Moreover, campaigns like 
PLANET/RoboNet intensify the sampling over the peak region, 
which significantly enhances the chance to detect and character- 
ize super-Earths as compared to a standard sampling scheme. On 
the other hand, our simulations overestimate the chance to detect 
OGLE 2005-BLG-390Lb-like in low magnification events, due 
to the increased anomaly durations when a giant source sweeps 
over a planetary caustic, whereas a main sequence source star 
would on the contrary produce a signal with higher amplitude 
but shorter duration. 

Despite its simplicity, our simulation however shows ro- 
bustly that in any case OGLE 2005-BLG-390Lb-like planets 
have a modest but not negligible probability for discovery, while 
Jovian-type planets should be found much more easily. With 
OGLE 2005-BLG-390Lb already b eing the third reporte d de- 
tection of a planet by microlensing dBeaulieu et alj[2006b after 
two Jovian planets, we therefore have the first observational evi- 
dence for the suggestion that sub-Neptunes are common around 
M-dwarfs and much more frequent than Jovian compa nions. The 
detection of another super- Earth bv lGould et all J2006h gives fur- 
ther support to this finding. 



of q = 7.6 x 10~ 5 is about 1%, whereas it becomes about 50% 
for q = 4.3 x 10 -3 , resembling a Jupiter-mass planet around the 
same host st ar, i.e. the sensit i vity to Jupiters in this case is, as was 
reported by iBeaulieu et al.1 (|2~006). about 50 times higher than 
it is to sub-Neptunes. Unsurprisingly the sensitivity to Jovians 
only shows a weak dependance on the applied detection criteria, 
since the majority of the associated anomalies are strong enough 
to be detected even by our most stringent detection demands. 
The sensitivity for OGLE 2005-BLG-390Lb-like mass ratio on 
the other hand is much more affected by the choice of the de- 
tection threshold and, when integrated over the lensing zone, 
ranges from ~ 0.5 — 3 % for criteria C3 to C\. For compari- 
son we have also computed and plotted the detection efficiency 
using a pure A\ 2 threshold with A\ 2 = 60, 220, 400, assum- 
ing the same light curve sampling rate (hourly) and photometric 
accuracy (1 %). We note that whereas the efficiencies for Jovian 
planets again remain rather unaffected at levels of ~ 50 — 60%, 
the detection efficiencies drops to 2—5 % for OGLE 2005-BLG- 
390Lb-like companions. In fact, a criterion assuming A\ 2 = 60 
is much less demanding that the C2 criterion, since the latter 
rejects a lot of anomalies which are too short lived to pass the 
minimal sampling requirement of 15 points. 

Averaged over the lensing zone one finds sensitivity ratios 
between Jupiters and sub-Neptunes of up to 25 for OGLE 2005- 
BLG-390-like time scales of t-g = 11.0 days. For time scales 
of around tE = 22.4 days, which is the median time scale of 



6. Summary and conclusions 

We have re-examined the photometric data on OGLE 2005- 
BLG-390Lb to look for traces of additional Jovian companions, 
finding that additional planets more massive than 1 Mj and in 
orbits of ~ 1.1 — 2.3 AU would have caused a detectable signal 
in more than 50% of the cases, which however was not observed. 
Planets with masses of 3 Mj and above between ~ 1.5 — 2.2 AU 
would have revealed themselves in the data with a probability of 
70%. 

Planet formation model s based on sequential accretion pro- 
cesses by llda & Linl d2005l) predict that the creation of gas gi- 
ants is strongly suppressed around M-dwarfs for practically the 
whole range of their model parameters. In agreement with such 
theoretical predictions we do not find in our data any indication 
of an additional companion to OGLE 2005-BLG-390Lb in the 
Jovian mass regime, however we are not able to definitely ex- 
clude this possibility. 

Assuming the natural unfiltered uniform distribution of lens- 
source impact parameters above an un-blended magnification 
threshold of A = 1.34 and a simple sampling pattern, we find a 
detection probability for an OGLE 2005-BLG-390Lb-like planet 
in events involving giant source stars of 1 — 3% at its angular sep- 
aration of 1.61 times the angular Einstein radius 6*e, whereas the 
average over the lensing zone, i.e. separations between 0.62 and 
1.62 9e, becomes 2 — 5%. With detection probabilities of a few 
percent, our discovery of OGLE 2005-BLG-390Lb provides the 
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first observational evidence in support of the prediction by cur- 
rent planet formation theories that sub-Neptune mass planets are 
common. 
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Appendix A: Conversion to physical lens 
parameters 

The parameters that can be extracted from the observed light 
curve are not sufficient for directly determining the properties of 



the planet and its host star, such as their masses, the orbital radius 
and period, as w ell as their distance. However, as discussed by 
iDominild d2006l) . one can derive probability densities by means 
of Bayes' theorem under the assumption of a mass spectrum for 
the lens stars, the lens and source distance, following the spa- 
tial mass density of stars in the Milky Way, and their velocity 
distribution. The planet itself is characterized by the planet-to- 
star mass ratio q = m/M and the separation parameter d, where 
d 9e is the instantaneous angular separation between the planet 
and its host star. All further information relies on the time-scales 

= $e//^ an d t+ = during which the source moves by 

the angular Einstein radius 6*e or its own angular radius 6>* rela- 
tive to the lens on the sky, where /i denotes its relative proper mo- 
tion. With 0* determined from the source magnitude and color, 
one obtains /i and 6*e- 

With Dl and Dg denoting the lens or source distance, re- 
spectively, let us define fractional distances x = D^/Dg and 
y = D S /R GC , where R GC = (7.62 ± 0.32) kpc is the dis- 
tance to the Galactic center dEisenhauer et al.l [2005b . If one ig- 
nores selection effects for the source stars related to their lu- 
minosity function as well as extinction, the probability of find- 
ing these between D$ and D$ + dD$ becomes proportional to 
Dg ps{Ds), where ps(Ds) denotes the volume mass density of 
the source stars, while pl(-Dl) denotes the volume mass density 
of the lenses. 

We further assume that the lens star mass spectrum 
$m/Mq (M/M® ) stretches from M m - ln and M max and does not 
depend on the lens distances, while we consider the distribu- 
tion of the effective lens velocity v = Z?l A 4 to depend on both 
the lens and source distance. With a characteristic velocity v c , 
the dimensionless velocity parameter £ = v/v c is distributed as 
^c(C,x,y) = v c <S> v (v c (,x,y). 

With 
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the measured values of and /i determine the detection effi- 
ciency as function of the projected separation r = d 6>e and 
the planet mass m as 

JVW/M 

e r ,m(r,m) = J J x 

mm(m/MQ,M m in/M®) 



y R GC m s ^e,0 \ M/M® J ' M/M® 
x PM/M Q (M/MQ,y;i 1 g E ,r] t ,)dyd(M/M®) , 



(A.4) 



where p M /Mq denotes the probability density of the mass of the 
primary, which reads 
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x I 2 i Trm I PsiyRcc) , (A.5) 



where 7V = // Pm/m q (M/Mq , y; ^ E , T]„) dy d(M/M Q ). 

The function f y (y) denotes a prior for the source distance, 
for which we adopted 

/tW = -t^= — ex p i — r ( A - 6 ) 



so that lny follows a normal distribution around In?/ ) with 
standard deviation <Ti n y = a y /y. 

With the assumption of circular orbits, de-projection assum- 
ing random orientation and phase, finally yields the detection 
efficiency as function of the orbital radius a and the planet mass 



(a,m) — / £f. m (y/ 1 — w 2 a, m) dw . (A.7) 



Fo r the mass and ve locity distribution of the Milky Way, we 
follow D ominikl (120061) in adopting a double-exponential disk 
and a barred bulge as well as mass spectra for disk and bulge 
lenses discussed bv lChabrierl (120031) . 

For OGLE-2005-BLG-390, we inferred i E = 11.0 d and 
t* = 0.28 d, so that with 9+ = 5.25 //as, one finds \i = 
18.6 /iasd" 1 and 6*e = 205 /ias. Typing of the source by 
iBeaulieu et all d2006h . yielded D { ° ] = 8.0 kpc and a Ds = 
2.0 kpc. 



